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ABSTRACT 

Ultra-luminous Compact X-ray Sources (ULXs) in nearby spiral galaxies and 
Galactic superluminal jet sources share the common spectral characteristic that 
they have unusually high disk temperatures which cannot be explained in the 
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framework of the standard optically thick accretion disk in the Schwarzschild 
metric. On the other hand, the standard accretion disk around the Kerr black 
hole might explain the observed high disk temperature, as the inner radius of 
the Kerr disk gets smaller and the disk temperature can be consequently higher. 
However, we point out that the observable Kerr disk spectra becomes significantly 
harder than Schwarzschild disk spectra only when the disk is highly inclined. 
This is because the emission from the innermost part of the accretion disk is 
Doppler-boosted for an edge-on Kerr disk, while hardly seen for a face-on disk. 
The Galactic superluminal jet sources are known to be highly inclined systems, 
thus their energy spectra may be explained with the standard Kerr disk with 
known black hole masses. For ULXs, on the other hand, the standard Kerr disk 
model seems implausible, since it is highly unlikely that their accretion disks are 
preferentially inclined, and, if edge-on Kerr disk model is applied, the black hole 
mass becomes unreasonably large (> 300M©). Instead, the slim disk (advection 
dominated optically thick disk) model is likely to explain the observed super- 
Eddington luminosities, hard energy spectra, and spectral variations of ULXs. 
We suggest that ULXs are accreting black holes with a few tens of solar mass, 
which is not unexpected from the standard stellar evolution scenario, and their 
X-ray emission is from the slim disk shining at super-Eddington luminosities. 

Subject headings: superluminal jet sources: ultra-luminous X-ray sources: accre- 
tion disks, slim disks: Schwarzschild black holes, Kerr black holes 


1. Introduction 

Ultra-luminous compact X-ray sources (ULXs) have been found in nearby spiral Galax- 
ies, with typical 0.5 - 10 keV luminosities 10 39 to 10 4 ° erg s -1 (e.g., Fabbiano 1988; Colbert 
and Mushotzky 1999; Makishima et al. 2000; Colbert and Ptak 2002; Foschini et al. 2002). 
These luminosities are too small for AGNs, and most ULXs are in fact located significantly 
far from the photometric center of the galaxies. On the other hand, ULXs are too luminous 
to be considered as the same class of the compact binary X-ray sources in our Galaxy whose 
luminosities are almost always < 10 39 erg s _1 . 

Significant time variations have been detected from ULXs, and their energy spectra 
are successfully modeled with emission from optically thick accretion disks (Okada et al. 
1998; Mizuno et al. 1999; Kotoku et al. 2000; Makishima et al. 2000; Mizuno, Kubota 
and Makishima 2001), as is the case for the “High” state (= Soft-state) of Galactic black 
hole candidates (e.g., Tanaka and Lewin 1995). In addition, discovery of bimodal spectral 
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transitions (Kubota et al. 2001) and orbital modulations (Bauer et al. 2001; Sugiho et al. 
2001) from several ULXs further demonstrate their resemblance with Galactic black hole 
candidates. These observational facts suggest that ULXs are moderately massive black holes, 
which may be scale-up versions of the Galactic black holes. So that the observed luminosities 
of ULXs, ~ 10 40 erg s -1 , do not exceed the Eddington limit Lem = 1.5 x 10 38 (M/M©) erg 
s _1 , the black hole mass has to be as large as or greater than ~100 M 0 , assuming isotropic 
emission. How to create such “intermediate” mass black holes, if truly exist, is an intriguing 
question (e.g., Ebisuzaki et al. 2001). On the other hand, if X-ray emission is unisotropic, the 
black hole mass is not required to be so large. For example, a beaming model for ULXs has 
been proposed (King et al. 2001; Kording, Falcke and Markoff 2002); however, presence of a 
bright nebula surrounding M81 X-9 (Wang 2002) suggests that the nebula is powered by the 
central X-ray source and that strong X-ray beaming is rather unlikely. Another possibility 
to solve the super-Eddington problem through unisotropic emission is a geometrically thick 
accretion disk (Watarai, Mizuno and Mineshige 2001), which we favorably consider later 
(section 4.1). 

GRS 1915+105 and GRO J1655-40 are the two well-known Galactic superluminal jet 
sources and established black hole binaries with reliable mass measurements (M « 7M© 
for GRO J1655-40 [Orosz and Bailyn 1997] and « 14M© for GRS 1915+40 [Greiner, Cuby 
and McCaughrean 2001]). If soft X-ray energy spectra of these sources are fitted with an 
optically thick accretion disk model, their characteristic disk temperatures are ~1.3 - 2.0 
keV (e.g., Belloni et al. 1997; Zhang et al. 1997; Zhang, Cui and Chen 1997; Tomsick et 
al. 1999). These values are systematically higher than those of ordinary and well-studied 
soft-state black hole candidates such as Cyg X-l and LMC X-3, whose disk temperatures 
are almost always less than ~ 1 keV (e.g., Tanaka and Lewin 1995). 

Okada et al. (1998) found that a luminous ULX in IC342 (~ 2 x 10 40 erg s -1 ; “Source 
1”) has unusually high disk temperature (~1.7 keV), which is similar to the accretion disk 
spectra of Galactic superluminal jet sources, and pointed out that such a high disk tem- 
perature cannot be explained in the framework of the “standard” accretion disk around a 
Schwarzschild black hole (see section 2.1 and 3.2). The “standard” disk denotes the situa- 
tion that all the gravitational energy release is converted to thermal radiation, in contrast 
to the “slim” disk which is an optically and geometrically thick disk with dominant energy 
advection (section 4.1). Makishima et al. (2000) summarized ASCA observations of seven 
ULXs, and concluded that the unusually high accretion disk temperature is a common spec- 
tral property of ULXs. Zhang, Cui and Chen (1997) and Makishima et al. (2000) suggested 
that the unusually high disk temperature of the Galactic superluminal jet sources and ULXs 
may be explained if they harbor fast rotating black holes, since inner edge of the accretion 
disk gets closer to the black hole in the Kerr geometry and the disk can be hotter. On the 
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other hand, Watarai et al. (2000) and Watarai, Mizuno and Mineshige (2001) proposed that, 
instead of the standard disk, the slim disk model may explain the X-ray energy spectra of 
Galactic superluminal jet sources and ULXs. 

In this paper, we focus on the “too-hot disk” problem of ULXs and Galactic superluminal 
jet sources. We apply standard accretion disk models in Newtonian, Schwarzschild, and Kerr 
cases, and discuss how the relativistic effects and the disk inclination affect the black hole 
mass and the mass accretion rates obtained from the model fitting. We see that observed 
hard spectra of Galactic superluminal jet sources may be explained by strong relativistic 
effects in the highly inclined Kerr accretion disk. On the other hand, the hard spectra 
and the super- Eddington problem of ULXs are difficult to explain in the framework of the 
standard accretion disk model. We see that the slim disk model may explain the observed 
super- Eddington luminosities and hard energy spectra of ULXs more naturally. 


2. Characteristics of the Standard Accretion Disk Model 

Before quantitative discussion of the observed energy spectra and accretion disk param- 
eters, we summarize important observational characteristics of the standard accretion disk 
model. 


2.1. Disk Temperature 

Let’s assume a standard optically thick accretion disk (Shakura and Sunyaev 1973) in 
the Schwarzschild metric, in which case the last stable orbit around the non-rotating black 
hole is 6 r g , where the gravitational radius r g = GM/c 2 is defined. We identify the last 
stable orbit with the innermost disk radius, r in . The radial dependence of the “effective” 
temperature of an optically accretion disk may be written as, 

m , . ( 3GMM _ , , ,1 1/4 

Teff(r) - | 8jrar 3 R R(. r / r in)j » ( 1 ) 

where Rr is ^1 — y/r in /r^ in the Newtonian case, otherwise includes additional general 
relativistic correction (see e.g., Krolik 1999). The effective temperature is zero at r = r in , 
peaks at r « 8 r g , and decreases with oc r~ 3 / 4 outward. 

Because of the Comptonization, the color temperature of the local emission becomes 
higher than the effective temperature, and their ratio is almost constant at « 1.7 (section 
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2.3). Therefore, the observable disk color temperature also peaks at r « 8r g , and takes the 
maximum values, 
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where we define the critical mass accretion rate Me = 2.9 x 1O 18 (M/M 0 ) g s _1 , so that 
M = Me gives the Eddington luminosity (see appendix A). The small temperature variation 
reflects the inclination effects (from face-on to i = 80°), such that inclined Schwarzschild disks 
have slightly harder spectra (section 2.4 and Appendix B). 

From equation (2), we can see that an optically thick accretion disk around a 7 M 0 or 
100 Mq Schwarzschild black hole may not have a higher color temperature than ~1.3 keV or 
~0.7 keV, respectively, unless super-Eddington mass accretion is allowed. To explain higher 
disk temperatures, either super-Eddington mass accretion rates or unreasonably small mass 
is required. This is the “too-hot disk” problem of the Galactic superluminal jet sources and 
ULXs we are concerned in this paper. 


2.2. Multi-color disk approximation 

We point out that the maximum disk color temperature as expressed with (2) is a 
quantity directly constrained from observations. If we approximate optically thick accretion 
disk spectra with a simple multicolor disk blackbody model (MCD model; Mitsuda et al. 
1994; Makishima et al. 1996) in which local emission is a blackbody and radial dependence 
of the temperature is simplified as T(r ) = T in (i^ n /r) 3 ^ 4 , the apparent inner disk radius 
Ri n and temperature T in will be the independent model parameters 2 , such that the spectral 
shape is determined only by Tj n , and the flux is proportional to R 2 n T^ n . When observed 
spectra are fitted with the MCD model, T in may be identified with the maximum disk color 
temperature (2), while the true inner disk radius rj n and apparent radius Ri n have a relation 
n n « 0.4 (Tcoi/Teff) 2 Ri n (Kubota et al. 1998). 


2 In this paper, rj„ denotes the real disk inner radius, while Ri n is the apparent inner radius as a MCD 
model parameter. Appendix C gives conversion formulae between the MCD parameters {R in and 7i„) and 
realistic disk parameters ( M and M) in the Schwarzschild case. 
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2.3. Local Energy Spectra 

In the inner part of the optically thick accretion disk, the electron scattering opacity is 
dominant, and the hot disk atmosphere distorts the emergent spectra via Comptonization. 
Such precise theoretical calculation of the accretion disk spectra for Galactic black hole 
candidates has been made by several authors (e.g., Shimura and Takahara 1995; Ross and 
Fabian 1996; Blaes et al. 2001). These authors basically agree that the local spectrum from 
each ring of the disk can be approximated by the “diluted black body”, ( T e ff/T coi ) 4 B(E , T co i), 
where B(E,T co i ) is the Planck function with color temperature T co i (> T e //), in the ~0.5 - 
10 keV band (where most observations are relevant). Also, ratio of the color temperature to 
the effective temperature, Tco(/^e//, is virtually constant at 1.7 - 1.9 along the disk radius. 
Hence, equation (2) is valid, and simple MCD approximation is useful to describe observed 
disk spectra. Throughout this paper, we use Tcoi/T e ff = 1.7 to fit the observed spectra. 

It should be noted that the mass obtained by fitting the standard accretion disk model 
to the observed spectra is proportional to ( T co i/T e ff ) 2 , as long as the diluted black body 
approximation is valid. This can be readily seen as follows: The accretion disk spectrum 
may be expressed as, 

r 2 * r (Iff) b ( e > m - dr 

&r in 2 [ 2?r(r/r in ) ( 7 ^] B (E,T col (r/r in , M, M)) d(r/r in ) 

J(r/r in )= 1 \lcolJ V ' 

= r°° 2w(r/r in )B(E,T e!l (r/r ln ,(T'„/Trt) 2 M,M)) dfr/r*,), 

V C / J(r/r in )= 1 V ' 

(3) 

where we have used Ti n = 6GM/c 2 and T coi oc M 1 / 4 (r/r in ) 3 7/4 ((T e ///T C0 ;) 2 M) from 
equation (1). Equation (3) indicates that the accretion disk spectrum with the black hole 
mass M and the diluted black body local spectrum is identical to the disk spectrum with 
the black hole mass ( T e ff/T co i ) 2 M and the black body local spectrum. 


2.4. Relativistic Effects 

Due to the relativistic effects in the Schwarzschild metric, observed disk spectral shape 
is mildly inclination-angle dependent, such that inclined disk spectra become slightly harder 
(appendix C, figure 1). This is due to the gravitational redshift and Doppler boosts near the 
inner edge of the accretion disk. Note that the relativistic effect is not significant outside 
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of 6 r g , so that the Schwarzschild disk spectra are not very much different from those of 
Newtonian disks. This is in contrast to the extreme Kerr case, in which the relativistic effect 
is enormous near rj n « r g and the disk spectrum is extremely inclination dependent (see 
below; figure 1 and 2). 

While in the Schwarzschild metric the last stable orbit around the black hole is 6 r p , in 
the Kerr metric it can go down to 1.24 r g with an extreme angular momentum of a = 0.998. 
As the inner edge of the accretion disk approaches the black hole, more gravitational en- 
ergy is released, and the inner disk temperature can get higher. X-ray energy spectra from 
optically thick Kerr accretion disks have been calculated by many authors including Cun- 
ningham (1975), Connors, Piran and Stark (1980), Asaoka (1989), Sun and Malkan (1989), 
Laor, Netzer and Piran (1990), Hubeny et al. (2000) and Gieliriski, Maciolek-Niedzwieck and 
Ebisawa (2001). In this paper, we calculate the Kerr disk spectral model using the transfer 
function by Laor, Netzer and Piran (1990) for a — 0.998. To make a comparison easier, 
we assume r in = 1.24 r g , and the same local spectrum as we did for the Schwarzschild disk 
model (section 2.3), namely diluted black body with T co i/T e ff = 1.7. 

In Figure 1, we show comparison of the Schwarzschild and Kerr disk spectra, as well as 
Newtonian one, for the face-on and a highly inclined case ( i = 80°). The biggest difference 
between the Kerr disk model and the other two models is the presence or absence of the 
innermost region of the disk from 1.24 r g to 6 r g . When the disk is close to face-on, the 
contribution from this part is not significant because of the gravitational red-shift and light- 
bending; as a result, the total Kerr disk spectrum is not very much different from the 
Newtonian or Schwarzschild ones. On the other hand, emission from the innermost part 
becomes very significant for a highly inclined disk due to Doppler boosting. Consequently, 
a near edge-on Kerr disk spectrum comprises a significant amount of high energy X-ray 
photons compared to the Newtonian or Schwarzschild disks. 

Inclination angle dependence of the observed Kerr disk spectrum is different in different 
energy bands (figure 2). At lower energies, where most of the emission is from the outer parts 
of the disk (r > 400 r g ), the flux is proportional to the projected disk area and decreases 
with inclination, just as Newtonian disks. On the other hand, higher energy fluxes from the 
innermost parts of the disk (1.26 r g < r < 7 r g ) are enhanced as the disk is more inclined, 
due to strong Doppler boosting. Consequently, inclined Kerr disks are brighter than the 
face-on disks in the highest energy bands, as previous authors have already pointed out. 


- 8 - 


3. Application of the Standard Disk Model 

We apply the standard accretion disk spectral models to ASCA archival data of GRO 
J1655-40 and IC342 Source 1. The observations were made on 1995 August for GRO J1655- 
40 (when the source is in a bright state) and 1993 September for IC342 Source 1. These 
datasets may be considered “exemplary”, so that the same spectral data of GRO J1655-40 
have been analyzed in Zhang et al. (1997), Zhang et al. (2000) and Gielihski, M., Maciolek- 
Niedzwiecki, A. and Ebisawa, K. (2001), and those of IC 342 Sourcel in Okada et al. (1998), 
Mizuno et al. (1999), Watarai, Mizuno and Mineshige (2001) and Kubota et al. (2001). 

We assume the diluted blackbody local spectra with T co i/T e ff = 1.7 in sections 3.1 and 
3.2. Difference among the Newtonian, Schwarzschild, and Kerr disk models are only the 
relativistic effects, including the variation of the innermost radius and energy conversion 
efficiency (appendix A). The Schwarzschild disk model we use (GRAD model) is explained 
in Hanawa (1989) and Ebisawa, Mitsuda and Hanawa (1991) (see also appendix B). The 
GRAD model is available in the XSPEC spectral fitting package, and this model also provides 
the Newtonian disk spectra, when the relativistic switch is turned off. 


3.1. Application to GRO J1655-40 

In Table 1, we show results of spectral fitting of the ASCA GIS spectra with Newtonian, 
Schwarzschild, and extreme Kerr disk models. System parameters have been determined 
from optical observations (Orosz and Bailyn 1997), such that distance is 3.2 kpc, inclination 
angle is 70°, and M = 7M Q . We fixed the distance, while several different inclination angles 
are tried from 0° to 80°. The energy spectrum has a power-law hard-tail but contribution of 
which is minor below 10 keV and hardly affects discussion of the accretion disk spectrum. 
The power-law slope is fixed at 2.5 which matches the simultaneous high energy observation 
with BATSE (Zhang et al. 1997), and its normalization is allowed to be free within the 
range acceptable by BATSE. Hence, the free parameters are M, M, power-law normalization, 
and the hydrogen column density. Besides minor local features (see Gielihski, Maciolek- 
Niedzwieck and Ebisawa 2001), all the models with different inclination angles can reasonably 
fit the observed accretion disk spectrum. 

Note that the Schwarzschild disk model at the correct inclination angle (i = 70°) gives 
a too small mass (1.8 M & ) which is unacceptable (red curve in Figure 3 left). If we fix the 
mass at the correct value (7 M©), the disk spectrum cannot be hard enough to explain the 
observed spectrum (the “best-fit” is shown with green curve in Figure 3 left): The “too- hot” 
disk problem is evident (equation 2). 


-9- 


It is of interest to see how the black hole mass depends on the inclination angle. Ob- 
servation can constrain the projected disk area, which is proportional to square of the mass. 
Hence M oc l/\/cosi in the Newtonian case, and the black hole mass at i = 80° is 2.4 times 
higher than that for the face-on disk. On the other hand, the mass increasing factor is 2.9 
in the Schwarzschild case, and as much as ~ 10 in the extreme Kerr case. This is because 
inclined relativistic disk spectra become harder, which is most conspicuous in the extreme 
Kerr disk (section 2.4; figure 1). Given the observed spectra, to compensate the spectral 
hardening with inclination, the mass becomes necessarily larger (equation 2). 

The extreme (a=0.998) Kerr disk model with the correct inclination angle gives M = 
15.9 M q . Compared to the fit with Schwarzschild disk model (1.8M 0 ), significant increase of 
the mass does indicate the spectral hardening of the Kerr disk model. In fact, the derived 
mass is too large to be consistent with the realistic mass 7 Mq, which probably suggests that 
a=0.998 is too high. In fact, Gielinski, Maciolek-Niedzwieck and Ebisawa (2001) applied the 
Kerr disk model to the same ASCA energy spectrum of GRO J1655-40, and concluded that 
a is likely to be between 0.68 and 0.88 to be consistent with M = 7 Mq. To summarize, in 
agreement with Gielinski, Maciolek-Niedzwieck and Ebisawa (2001), a standard Kerr disk 
model may explain the observed accretion disk spectra of GRO J1655-40 at the inclination 
angle i = 70°. The Kerr metric is required, but not with the maximum angular momentum. 

Probability that GRO J1655-40 has a spinning black hole is also suggested by the recent 
discovery of 450 Hz QPO (Strohmayer 2001), as this frequency is higher than the Kepler 
frequency at the innermost stable orbit of a 7 M 0 Schwarzschild black hole. Precise analysis 
of the QPO characteristics from GRO J1655-40 suggests that the black hole is in neither 
a Schwarzschild nor a maximal Kerr black hole (Abramowicz and Kluzniak 2001), which is 
consistent with our energy spectral analysis. 


3.2. Application to IC342 Source 1 

We assume the distance to the source 4 Mpc (Okada et al. 1998 and references there 
in). If isotropic emission is assumed, the luminosity will be 1.7 x 10 40 erg s -1 (1 - 10 keV) 
at this distance, thus M > 100M e is expected so as not to exceed the Eddington luminosity. 
In table 2, we summarize results of the spectral fitting of the ASCA GIS spectrum with 
Newtonian, Schwarzschild, and Kerr disk models, for different inclination angles. All the fits 
are acceptable, thus models may not be discriminated based on the fitness for the observation. 
Note that the disk luminosities always exceed the Eddington limit more than 10 times in the 
Schwarzschild case. 
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The right panel in figure 3 indicates the ASCA GIS spectrum and the best-fit Schwarzschild 
disk model (in red). For comparison, we show a Schwarzschild disk spectrum with M = 
100 Mq at the Eddington limit (in green). We can clearly see that such an accretion disk has 
a too low temperature to explain the observed hard spectrum. 

The face-on Kerr disk model gives M = 27.3 Mq and M = 2.0 x 10 20 g s -1 (M = 16 Me). 

A factor of ~ 3 increase of the mass compared to the Schwarzschild case is due to slight 
hardening of the face-on Kerr disk spectrum. Note that the super-Eddington problem does not 
disappear, as the face-on Kerr disk spectrum is not very different from the Schwarzschild one 
(section 2.4; figure 1). If we assume a very inclined disk with i = 80°, we obtain M = 332 M 0 
and M = 1.5 x 10 20 g s -1 (M = 1.0 Me)- Now the super-Eddington problem is solved, that 
is consequence of the fact that the inclined Kerr disk spectrum is much harder than the 
Schwarzschild one. However, there will be two serious problems to accept the inclined Kerr 
disk model for ULXs in general: First, we do not know a mechanism to create ~ 350 M 0 
black holes. Second, it is very unlikely that most of the accretion disks in ULXs are largely 
inclined when seen from the earth. 


3.3. Examination of Variants of Standard Disk Model 

In this section, we critically examine two possibilities which have been proposed to make 
the standard accretion disk spectra look harder. 


3.3.1. Large T^i /T e ff? 

A naive solution to explain the apparently hard disk spectrum may be to allow T co i/T e ff 
to be much greater than the standard value ~ 1.7, in which case the color temperature of 
the disk can be much higher for the same mass and mass accretion rate, while the effective 
temperature remains the same (equation 2). In other words, for a given observed disk flux 
and spectrum, M a (T coi /T e //) 2 , thus a larger mass is allowed (equation 3). 

In fact, in the case of GRO J1655-40, (Section 3.1, table 1), if T co i/T e ff = 1.7-\/7Mq/1.8Mq 
3.4, the observed spectrum is explained with the Schwarzschild disk model with 7 M© for the 
same distance and inclination angle. Borozdin et al. (1998) reached a similar conclusion that 
Tcoi/Teff — 2.6 is required for GRO J1655-40 to fit with a Newtonian disk model whose 
inner disk radius is 3 times the Schwarzschild radius. For IC342 Source 1, we need a still 
higher value of T co i/T e ff = U1/IOOM0/8.9M© =5.7, to fit with a face-on Schwarzschild 
disk model around a 100 Mq black hole. 
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Although having high T co i/T e ff values might seem to solve the “too-hot” disk problem, 
T co i/T e ff ^ 2 for Galactic sources and ULXs is most certainly excluded both from theory and 
observation. Independent theoretical calculations agree with values of = 1.7 - 1.9 

for accretion disks around ~ 10M©black hole (Shimura and Takahara 1995; Ross and Fabian 
1996; Blaes et al. 2001). Values of T col /T e ff increase very slowly with mass, but cannot be 
higher than ~ 2.5 even for a 10 6 M© black hole (Ross, Fabian and Mineshige 1992; Blaes et 
al. 2001). From observations, when applied to accretion disk spectra of Galactic black hole 
candidates or weakly magnetized neutron stars, T co i/T e ff = 1.7 -1.9 gives reasonable mass 
values which are consistent with those derived from more reliable methods (e.g., Ebisawa, 
Mitsuda and Hanawa 1991; Ebisawa et al. 1993; Shimura and Takahara 1995; Dotani et al. 
1997). If we happen to use T C0 (/T e // « 2.6 for an accretion disk around a weakly magnetized 
neutron star, the neutron star mass exceeds ~ 3M©, which is not acceptable. 


3.3.2. Comptonization of disk photon? 

If a standard optically thick disk is shrouded by a thermal corona whose temperature 
is slightly higher than the disk temperature, significant parts of the disk photons are comp- 
tonized and appear in higher energy bands. Thus, in such a situation, the observed hard 
energy spectra from Galactic superluminal jet sources and ULXs might be explained. 

The comptonization model may phenomenologically fit the observed spectra of GRO 
J1655-40 (Zhang et al. 2000; Zycki, Done and Smith 2001). Fitting with a thermal comp- 
tonization model (in which seed photons have a single temperature black body spectrum) 
gives typically the black body temperate < 0.5 keV, plasma temperature ~ 1.5 keV and 
scattering optical depth ~10. At such a low temperature and high optical depth, the heavy 
elements may not be fully ionized, and the plasma can be optically thick to photoelectric 
absorption. Presumably, such a physical situation is more reasonably described by standard 
optically thick accretion disk in which emergent spectrum is distorted by Comptonization 
(T co i > T e ff), which successfully explain the observation (section 3.1). 

We applied the comptonization model to IC342 Source 1. When the seed photon is 
assumed to be disk blackbody, we obtain T, n « 1.2 keV, plasma temperature ~ 2.5 keV, and 
scattering optical depth ~ 10 (x 2 /dof=51.8/46). Note that the T in cannot go down to ~ 0.7 
keV as desired for the Schwarzschild disk around a 100 M© black hole. Therefore, the original 
“too-hot” disk problem may not be solved even introducing the ad hoc comptonization 
plasma to harden the standard disk spectrum. When the seed photon spectrum is assumed 
to be black body, the black body temperature is 0.52 keV, the plasma temperature is 1.8 
keV and the optical depth is ~ 18 (x 2 /dof=51.8/46). Although this model may fit the data, 
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it is hard to interpret these parameters in the physical context, and the super- Eddington 
problem is unanswered. Presumably, the slim disk model explained in the next section is 
more likely for ULXs. 


4. The Slim Disk Model 


4.1. Characteristics of the slim disk model 


So far, we have considered the standard optically thick accretion disk (Shakura and 
Sunyaev 1973) in which radial energy advection is neglected and all the gravitational en- 
ergy release is converted to thermal radiation. In accretion disk theory (for a review, e.g., 
Kato, Fukue and Mineshige 1998), there is another stable optically thick solution, which 
is called “optically thick ADAF (advection-dominated accretion flow)” disk or, “slim” disk 
(Abramowicz et al. 1988), which takes place when mass accretion rate is extremely high. 
In the slim disk, all the gravitational energy release is not converted to the thermal radia- 
tion, but significant part of the energy is carried inward due to radial advection. Slim disk 
is geometrically thick, and can be much hotter than the standard disk (Kato, Fukue and 
Mineshige 1998; Watarai et al. 2000). 

One of the most significant observational characteristics of the slim disk is that the disk 
luminosity can exceed the Eddington limit by up to ~ 10 times (Abramowicz et al. 1988; 
Szuszkiewicz, Malkan and Abramowicz 1996; Kato, Fukue and Mineshige 1998). This may 
be qualitatively understood as follows: At any disk radius, local radiation pressure may not 
exceed the vertical gravitational force, such that 


F(r)< 


cGMh 

9 J 

kt z r 


( 4 ) 


where F(r ) is the energy flux, r the disk radius and h the half-thickness. Therefore, 


Ldisk — 


pf'out 
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k 




~ Lsdd 
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where r in and r out are inner and outer disk radius, respectively, such that In (r out /r in ) fa 10. 
In the standard disk, which is geometrically thin, h/r < 0.1, thereby Ldisk < ^Edd- On the 
other hand, with the slim disk in which h/r 1) Ldisk <10 Lsdd- 
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In addition to that super-Eddington luminosities are permissible, slim disk has following 
observational characteristics (e.g., Kato, Fukue and Mineshige 1998; Watarai et al. 2000; 
Watarai, Mizuno and Mineshige 2001): (1) As the mass accretion rate exceeds the critical 
rate Me, energy conversion efficiency decreases due to advection (appendix A, figure 7). 
Hence the disk luminosity is no longer proportional to the mass accretion rate, but saturates 
at ~10 Lsdd- (2) Innermost disk radius can be closer to the black hole than the last stable 
orbit, as mass accretion rate increases. Even in the Schwarzschild metric, innermost disk 
region inside 6 r g can emit significant amount of X-rays. Thus slim disk can be “hotter” 
than the standard disk. (3) If radial dependence of the disk effective temperature is written 
as T e ff oc r~ p (“p-free disk” below), the exponent p reduces from 0.75 (which is expected 
for the standard disk) to 0.5 as advection progressively dominates. Thus spectral difference 
from the standard disk will be noticeable. 


4.2. Slim disk for Galactic Sources? 

Watarai et al. (2000) suggested that the slim disk model may explain the apparently hard 
energy spectra of Galactic superluminal jet sources. In fact, energy spectra of GRO J1655-40 
(Kubota 2001), GRS 1915+105 (Yamaoka 2001) and XTE J1550-56 (Kubota 2001) are not 
fitted with a standard accretion disk model when their disk luminosities and temperatures 
reached maxima, but better fitted with the p-free disk model with p between 0.75 and 0.5. 
This is considered to be an evidence of emergence of the slim disk, when mass accretion rates 
are extremely high in these sources. 

On the other hand, when the disk luminosity and temperature are lower, GRO J1655- 
40 is considered to embrace the standard disk, since the standard disk model can fit the 
observed spectra well, and the inner disk radius is fairly constant over a large luminosity 
variation (Sobczak et al. 1999; Kubota, Makishima and Ebisawa 2001). Without advection, 
the standard Kerr disk with high inclination is successful to explain the observed hard disk 
spectra (section 3.1). Therefore, except when the disk luminosity is near the peak and the 
slim disk presumably takes place, the standard disk is expected to be present in the Galactic 
superluminal jet sources. 


4.3. Application to IC342 Source 1 

As for ULXs, the slim disk model seems to be more likely, as ULXs are intrinsically bright 
systems. Mizuno, Kubota and Makishima (2001) studied spectral variations of several ULXs, 
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and commonly found anti-correlation between the MCD parameters Ri n and T in , which is 
another observational feature of the slim disk (Watarai et al. 2000). Watarai, Mizuno and 
Mineshige (2001) calculated energy spectra of slim disks around intermediate mass black 
holes, and suggested that ULXs are slim disks around 10 - 30 M© black holes shining more 
luminously than Eddington limits. 

In the present paper, we fit the observed energy spectra with the same slim disk model 
calculated by Watarai, Mizuno and Mineshige (2001). Namely, a = 0.01, and the local emis- 
sion is a diluted blackbody with T co i/T e ff = 1.7. The model assumes the face-on geometry, 
and the distance is 4 Mpc. The energy spectra are calculated for each grid point of M = 10, 
32, 100 M© and M/(T£dd/c 2 ) = 1) 3, 10, 32, 100, 320, 1000. We fitted the observed spectra 
using XSPEC spectral fitting package with this grid-model, so that XSPEC interpolates the 
model spectra for other M and M values. 

Fitting result for IC342 Source 1 ASCA GIS spectrum is shown in table 3. We find 
M = 19.8 Mq and M/Mc ~ 20. This mass accretion rate would give 20 times the Eddington 
luminosity in the case of the standard disk. However, since slim disk is less efficient at 
M > Me (appendix A and figure 7), the disk luminosity is in fact only ~ 6 times the 
Eddington luminosity, which is allowed in the slim disk (section 4.2). 

The slim disk model can fit the observed spectra reasonably well (x 2 /dof = 1.40), but 
not as good as the standard disk models (x 2 /dof « 0.90; table 2) (figure 4). Also, the 
hydrogen column density is significantly larger with the slim disk model (~9x 10 21 cm -2 ) 
compared to those with the standard disk models (3-5 xlO 21 cm -2 ). This is because the 
spectral difference between the best-fit slim disk model and standard disk spectra is largest 
in the lower energy band (figure 5). We simply assumed diluted blackbody local spectra 
with constant T co i/T e fj , which is known to be a good approximation for the standard disk 
(section 2.3). However, precise theoretical calculation of the local spectrum in the slim does 
not exist yet, and more realistic slim disk spectral model is expected to fit the observed 
spectra better. 


4.4. Spectral Variation of IC342 Source 1 

The characteristic anti-correlation between Ri n and T jn discovered by Mizuno, Kubota 
and Makishima (2001) is considered to be a characteristic of the slim disk. Watarai, Mizuno 
and Mineshige (2001), through indirect comparison of the MCD parameters and slim disk 
spectra, claimed that this spectral variation may be explained with a slim disk a constant 
M only by varying the mass accretion rate. We demonstrate this claim by directly fitting 
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the ASCA IC 342 Source 1 spectra with the slim disk model. 

We use the same spectral datasets used by Mizuno, Kubota and Makishima (2001). The 
observation period is split into five periods depending on the flux levels. For each period, 
GIS and SIS spectra are fitted simultaneously to achieve the better statistics. In the top 
panel of figure 6, we show M and M variation obtained with the face-on GRAD model. 
We fixed the hydrogen column density at the average value (Ng = 5.2 x 10 21 cm -2 ), so 
that the free parameters are only M and M. We see clear anti-correlation between these 
two parameters, which is just rephrasing the Ri n -Ti n anti-correlation discovered by Mizuno, 
Kubota and Makishima (2001). Since M should not vary in reality, this relation is telling 
that the standard accretion disk is not a proper model for IC 342 Source 1. The spectral 
variation is clearly seen using only the 2-10 keV band, though less clear than using the 
entire 0.5 - 10 keV band. 

Since the present slim disk model does not fit the observed spectra well below ~1.5 keV, 
we study spectral variation only using 2-10 keV. For fair comparison, the hydrogen column 
density was fixed to the same average value as used in the GRAD model fitting. The M 
and M relation obtained from the slim disk fitting is shown in the bottom panel in figure 
6. We see that M in the range of 22.5 M Q to 23.9 M 0 is consistent with all the five spectra. 
Therefore, on the slim disk model, we may interpret the observed spectral variation of IC 
342 Source 1 as a result of simple mass accretion rate changes. 


4.5. Future Problems of the Slim Disk Model 

Compared to the standard disk model, X-ray spectral study from the slim disk has 
only a short history. As already stated above, the present slim disk model may not fit the 
observed ASCA spectra perfectly well, though the spectral variation is correctly described. 
However, there are several problems to calculate the slim disk spectra precisely, and this 
is not an easy problem. Although the face-on geometry is assumed in the present slim 
disk model, observed flux from a slim disk is considered to be strongly inclination angle 
dependent. The flux drops with inclination angle i more rapidly than cos i, because the 
horizontal photon diffusion time-scale gets longer than the in-fall time-scale (Kato, Fukue 
and Mineshige 1998). Also, if the slim disk exists around a Kerr black hole, the disk can be 
much hotter than the Schwarzschild case (Beloborodov 1998). In addition, effects of radial 
photon-trapping may not be negligible under extremely high accretion rates (Ohsuga et al. 

2002) , by which T coi /T e ff may have a significant radial dependence (Shimura and Manmoto 

2003) . More precise spectral model calculation from the slim disk is anticipated overcoming 
these difficulties. 
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5. Origin of the Ultra-Luminous X-ray Sources 

We have shown that the slim disk model is likely to explain the observed super-Eddington 
luminosities, hard energy spectra, and spectral variations of IC342 Source 1. Prom the spec- 
tral fitting, we obtained the black hole mass ~ 20 M 0 , and the disk luminosity ~ 6 times the 
Eddington luminosity. Thus, “intermediate mass black holes” are not required to explain 
the observed luminosity ~ 10 40 erg s -1 . 

In our Galaxy, the most massive stellar black hole is ~ 14 M e in GRS 1915+40 (Greiner, 
Cuby and McCaughrean 2001), as far as currently measured. However, more massive stellar 
black holes may well exist in other galaxies, since stellar evolution theory says main-sequence 
stars can have a maximum mass ~ 6OM 0 (Schwarzschild and Harm 1959). Grimm, Gilfanov 
and Sunyaev (2003) claims that there exists a universal luminosity function of the high mass 
X-ray binaries, such that the normalization is determined by the star-forming rate. Galaxies 
which are active in star formation tend to have more luminous and massive compact objects, 
and the maximum cut-off luminosity is ~ 10 40 erg s -1 . It is likely that ULXs with luminosities 
of ~ 10 40 erg s _1 are stellar mass black holes having a few tens of solar mass, which reside 
at the brightest end of the luminosity function. 


6. Conclusion 

1. In order to solve the “too-hot disk” problem of ULXs and Galactic superluminal jet 
sources, we have carefully studied energy spectra of standard accretion disk and slim 
disk. In particular, we have calculated the extreme Kerr disk spectral model, and stud- 
ied how relativistic effects and disk inclination angle affect the observed disk spectra. 

2. We have found that the Kerr disk model can successfully explain the observe hard 
spectra of GRO J1655-40, because the Kerr disk spectra become significantly harder 
than Schwarzschild ones when the disk is significantly inclined, that is the case for 
GRO J1655-40. 

3. The Kerr disk spectra are not significantly harder than the Schwarzschild disk when 
disk inclination angle is not large. We conclude that the standard Kerr disk model is 
not appropriate to explain the observed hard spectra of ULXs, as it is unlikely that 
accretion disks in most ULXs are preferentially inclined, and, if the inclined Kerr disk 
is applied, unreasonably large black hole mass (~ 300 M Q ) is required. 

4. We have calculated the slim disk spectra, and applied to the observed spectra of IC342 
Source 1. Although the current primitive slim disk model does not perfectly fit the 
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data, we found the slim disk can explain the observed super- Eddington luminosity, 
hard X-ray spectra, and spectral variation. In particular, the characteristic spectral 
variation is explained with constant mass at M ~ 20 M Q with simple mass accretion 
rate variation. 

5. We propose ULXs are a few tens of solar mass black holes, which reside at the bright 
end of the X-ray binary luminosity function. 
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XSPEC package. This research has made use of public data and software obtained from 
the High Energy Astrophysics Science Archive Research Center (HEASARC), provided by 
NASA’s Goddard Space Flight Center. This work was supported in part by Polish KBN 
grant 2P03D01718. 


A. Energy Conversion Efficiency and Critical Mass Accretion Rates 


Energy conversion efficiency of optically thick accretion disk ( 77 ) is dependent on disk 
model assumptions, so that a particular care is needed to interpret the disk parameters we 
obtain from spectral model fitting. In particular, conversion efficiency of the slim disk model 
is dependent on the mass and mass accretion rates (Watarai et al. 2000; Watarai, Mizuno 
and Mineshige 2001). 


In the Newtonian case, energy loss in the standard disk per second (= disk luminosity, 
Ldisk ) may be written as 

(GM L 2 


Ldisk — Af 


V n 


9 r - 2 

271 


GMM 
2 r in ’ 


where r, n is the innermost radius, and L is the specific angular momentum at r in , ^GMr in . 
If we take rj n as the last stable orbit in the Schwarzschild metric, 6 r g = 6 GM/c 2 , L disk = 
Ajlfc 2 , hence rj = 1/12 for the Newtonian disk. In the Pseudo-Newtonian potential, which 
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the present slim disk model adopts, 

, _ ( gm l 2 \ 

iV ‘ k ' r„, - 2r,, 2 r in 2 )’ 

where L = — 2 r p ) 2 . With r in = 6r g , we obtain 77=1/16. This is close enough 

to the correct efficiency in the Schwarzschild metric, 77=0.057. For the extreme Kerr case 
with a = 0.998 and r in = 1.24 r g , r) = 0.366. 

We may define the critical accretion rate Me so that Lsdd = ^Mcc 2 , where Lem Is the 
Eddington luminosity. For the Schwarzschild case with rj n = 6 r g> Me = 2.9 x 10 18 (M/M©) 
g s -1 , and for the Newtonian case with the same innermost radius, Me = 2.0 x 1O 18 (M/M 0 ) 
g s -1 . For the extreme Kerr case Me — 4.6 x 10 17 (M/M©) g s -1 . 

In the pseudo-Newtonian potential, if the innermost radius is constant at r in = 6 r g , 
Me = 2.7 x 10 18 (M/M©) g s 1 . However, in the present slim disk model, the inner radius 
becomes smaller than 6 r g , and the energy advection is dominant at high accretion rates 
(Watarai et al. 2000). Hence, the conversion efficiency is dependent on the mass and mass 
accretion rate. In figure 7, we show conversion efficiency of the slim disk for several mass 
and mass accretion rate values. We see that at low mass accretion rate limits, 77 « 1/16, but 
the efficiency significantly decreases with mass accretion rates when M > Me due to the 
advection effect. 


B. Bugs in the original GRAD model 

Three mistakes have been found in the appendix of Ebisawa, Mitsuda and Hanawa 
(1991; EMH) besides obvious typos: 


• (A8) in EMH should be read as follows: 

d9' ph _ cos i sin 0 P h + sin i cos 9 P h sin ip p h 
d9 ph ” sin 9' ph 

• (A15) in EMH should be read as follows: 

I- . _ L _ - 0 ' 5 [GMr%* 

£ L o V rc 2 \ rc 2 ) V rc 2 c dt ph 


(Bl) 


(B2) 


The bold-face part was 1 in EMH, which is wrong. 
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• (A22) in EMH should be read as follows: 


r col = 1.11 


Tcol/Te fA 

1.5 ) 


9irhj) 

0.15 




(B3) 


The bold-face part was 1.4 in EMH, which is incorrect. 


Below, equations (Al) to (A4) refer to those in the present appendix, not in EMH. 
Equation (Bl) was already fixed in the GRAD model code used in EMH, which was dis- 
tributed with XSPEC version ll.O.lae and previous ones. Equations (B2) and (B3) were not 
fixed in the code, and the wrong values were used in EMH. In addition, there was another 
bug in the code used in EMH: In the part to calculate the disk flux, 

F = 8.9038633D — 3 * F/Ratio * *4/(D/10.0) * *2 * (Em/1.4) * *2 
was wrong, and this has to be 

F = 8.9038633D — 3 * F /Ratio * *4/ (D/10.0) * *2 * Em * *2, (B4) 

where Em denotes the compact object mass. These bugs in equation (B3) and (B4) compen- 
sate each other (but not perfectly), which probably explains why we could not notice these 
bugs much earlier. The bug in equation (B2) affects only when the disk is inclined, and this 
effect is not significant. 

Besides the bug in equation (B2), one can see from equations (B3) and (B4) that the 
old GRAD spectra with M = 1.4 M 0 and the new GRAD spectra with M = 1.0 M e are 
identical. Therefore, besides the minor difference of redshifts which originates in (B2), when 
an observed spectrum was fitted with the old GRAD model, the mass of the compact object 
was erroneously 1-4 times higher than the correct mass given by the fixed GRAD model. This 
should be taken into account to interpret published results obtained using the old GRAD 
model (e.g., Ebisawa, Mitsuda and Hanawa 1991; Kitamoto, Tsunemi and Roussel-Dupre 
1992; Ebisawa et al. 1993; Dotani et al. 1997; Kubota et al. 1998). The mass accretion 
rate represented in the dimension of mass per time is unchanged. The fixed GRAD code is 
distributed with XSPEC version 11.1.0 and on. 

Note that the bug in the old GRAD code may compensate with an ambiguity of T co i/T e ff 
values. In fact, in Ebisawa, Mitsuda and Hanawa (1991) and Ebisawa et al. (1993), we used 
Tcoi/T e ff = 1.5 and obtained reasonable masses of black holes and neutron stars, which 
should be reduced by 1.4 times using the correct GRAD model. On the other hand, the 
mass estimate is proportional to ( T col /T e ff ) 2 (when distance is known) or (T col /T e fj) A (when 
d/y/M is known), thus using T co i/T e ff = 1.8 or 1.6 instead of 1.5 retrieves the original mass 
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values. The higher T c0 ;/T e // values are in fact consistent with a precise calculation carried 
out more recently (Shimura and Takahara 1995; Ross and Fabian 1996). 

We confirmed that the fixed GRAD model gives almost identical spectra with the 
Schwarzschild disk spectra calculated by Bhattacharyya, Bhattacharya and Thampan (2001). 
Only a minor difference is found when the accretion disk is more inclined than ~ 60°, when 
the model by Bhattacharyya, Bhattacharya and Thampan (2001) gives as much as ~ 10 % 
higher flux at ~ 10 keV. The discrepancy is probably due to the fact that the GRAD model 
does not take account of the photons whose trajectories change directions more than 180 
degrees. In fact, those photons may not be observable being blocked by the putative plasma 
near the black hole which is responsible for generation of the power-law hard-tail component. 


C. Comparison of Newtonian and Schwarzschild Accretion Disk Spectra 

Comparison of the GRAD and MOD model parameters have been made by EMH, which 
is affected by the bug explained above. Here, we present correct results using the new GRAD 
code after the bug is fixed. Also, we discuss how the relativistic effects (in the Schwarzschild 
case) affect the mass and mass accretion rate determination from observed accretion disk 
spectra. 

MCD model does not take into account the inner boundary condition, but simply as- 
sumes T(r) = T in {Rin/f) 3 ^, where R, n and T in are the apparent inner disk radius and tem- 
perature, respectively. MCD model is useful to represent observed spectra, since it has the 
two independent parameters which are directly constrained from observation, such that the 
spectral shape is determined only by T, n , and the normalization is proportional to Rf n cos i, 
where i is the inclination. Often disk inclination is unknown, so that Ri n is uncertain and 
proportional to (cos?) ~ 0 - 5 . 

Relationship between the MCD T in and R, n parameters and the mass and mass accretion 
rate in the Newtonian disk model with correct boundary condition is given as (see also 
Kubota et al. 1998), 



The multicolor disk spectrum and the Newtonian disk with correct boundary condition 
have almost identical spectral shape in 0.5 - 15 keV when equations (Cl) and (C2) are 
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considered. 


Formulae to relate the new GRAD model parameters and MCD parameters are the 
following: 
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(C3) 


(C4) 


The numerical coefficients have been obtained by fitting the new GRAD spectra with mul- 
ticolor disk model (see below) 3 . 


In Figure 8, we illustrate the relationship of the MCD model parameters and M and 
M in the Newtonian and GRAD models. These figures are made by simulating GRAD and 
Newtonian model spectra with various values of mass, mass accretion rate and inclination, 
using the Ginga response (sensitive in 1 - 30 keV), assuming the distance 1 kpc and exposure 
time 3000 sec exposure. These simulated spectra are fitted with the MCD model, and the 
best-fit Rin and T in are obtained. 

From equations (Cl) to (C4) and Figure 8, for a given observed spectrum for which Ri n 
and Ti n are determined, we can see the following (see also Table 2 and 1): 


• When fitting the same observed spectra, we will get a smaller mass and higher accretion 
rate with the GRAD model than with the Newtonian disk. 

• Differences in the mass and mass accretion rate between the GRAD model and the 
Newtonian models are the largest for the face-on disk, and get smaller with inclination. 


These characteristics of the GRAD model are due to relativistic effects, and may be 
understood intuitively as follows. Because the conversion efficiency of the GRAD model 
(0.057) is smaller than that of the Newtonian disk (0.083), the GRAD model requires a 


3 Equations (3) in Ebisawa, Mitsuda and Hanawa (1991) was wrong such that the old formula gives 1.4 
times larger mass for the same R ln and T in . 
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higher mass accretion rate to shine with the same luminosity. Thus we need a higher mass 
accretion rate with the GRAD model when fitting the same observed spectra. The observed 
disk flux is proportional to Lduk cos i, thus the disk luminosity induced from the observed flux 
is Ldisk oc (cos*) -1 . Consequently, the mass accretion rate is proportional to (cosi) -1 in the 
Newtonian case. In the GRAD model, this dependence becomes milder, because the face-on 
disk tends to be dimmer due to the light bending effect in the vicinity of the black hole, while 
those photons removed from the near face-on disks turn out to be the flux enhancement for 
the inclined disks. In fact, for the face-on disk, the GRAD mass accretion rate is ~ 1.8 time 
larger than that of the Newtonian disk, whereas it is ~ 1.1 times for i = 80°. 

Because of the gravitational redshift, the GRAD model gives a smaller characteristic 
disk temperature than the Newtonian disk with the same mass and mass accretion rate 
(compare equations C2 and C4). To achieve the same characteristic disk temperature of 
the Newtonian model, the ~1.8 times mass accretion rate increase is not sufficient, and 
the mass has to get smaller too. Since we measure the projected disk area oc M 2 cosi, 
the mass increases with inclination as (cosi) -0,5 in the Newtonian case, while the GRAD 
disk temperature gets higher because of the Doppler boost. Therefore, ratio of the GRAD 
mass to the Newtonian mass is the smallest for the face-on disk (~0.6) and gets larger with 
inclination (~0.73 for i > 60°). 

It is interesting that difference of the Newtonian and GRAD model parameters are the 
largest for the face-on disk and the smallest for highly inclined disks where relativistic effects 
are supposed to be most significant. This is because the gravitational redshift and Doppler 
boosts more or less cancel each other for a highly inclined Schwarzschild disk. This is in 
contrast to the extreme Kerr disk, where near edge-on disks show extremely hard spectra 
due to enormous Doppler boosts (see figure 1 and 2). 
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Table 1. GRO J1655-40 fit with Newtonian, Schwarzschild and Kerr disk models 
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Table 2. IC342 Source 1 fit with Newtonian, Schwarzschild and Kerr disk models 4 
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Table 3. IC342 Source 1 fit with Slim disk model 1 
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Comparison of Newtonian, Schwarzschild and Kerr disk spectra 



Fig. 1. — Comparison of Newtonian (black), Schwarzschild (red) and extreme Kerr (a = 
0.998; green) disk spectra for the face-on ( i = 0°) and a near edge-on disk (i = 80°). Inner 
disk radius is 6r g for the Newtonian and Schwarzschild disks, and 1.24 r g for the Kerr disk. 
The same mass accretion rate is assumed, which is so chosen to give the Eddington luminosity 
for the Schwarzschild disk. Other disk parameters are indicated in the figure. Note that the 
total disk luminosities are different for the three disk models because the energy conversion 
efficiencies are different (0.083, 0.057 and 0.366 for Newtonian, Schwarzschild and Kerr disks, 
respectively). 
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Kerr disk at d = 1 kpc, M= 1 M Q , Eddington limit 
fi = 0.9, 0.7, 0.5, 0.3 and 0.1 



o 0.01 0.1 1 10 


Energy [keV] 

Fig. 2. — The Kerr accretion disk spectra with an extreme angular momentum (a = 0.998), 
observed at the inclination angle /x = cos* = 0.9 (green; near face-on), 0.7 (yellow), 0.5 
(cyan), 0.3 (red) and 0.1 (black; near edge-on). Note the unit of the ordinate (keV 2 s" 1 
keV -1 cm -2 ) which facilitates to see the energy release per logarithmic energy. Solid lines 
indicate the total disk spectra, and contributions from inner (1.26r 9 < r < 7 r g ), middle 
(7 r g < r < 400 r g ), and outer parts ( 400 r g < r ) are plotted separately either by dotted line 
or broken line. The distance and mass are assumed to be 1 kpc and 1 M© respectively. The 
Eddington luminosity is assumed, and T co i/T e ff — 1. 




GRO J1655-40 


IC342 Source 1 



Fig. 3. — ASCA GIS energy spectra of GRO J1655-40 (left panel) and IC342 Source 1 
(right). The best-fit Schwarzschild disk models are indicated in red, in which case the best- 
fit masses are too small (1.8 M© and 8.9 M© for GRO J1655-40 and IC342 respectively) 
compared to the mass determined from optical observations (7 M© for GRO J1 655-40), or 
that expected from the observed luminosity (100 M© for IC342). Schwarzschild disk spectra 
with the expected masses give too low temperatures to explain the observed accretion disk 
spectra for both cases (shown in green). Additional power-law component (blue in left panel) 
is required to fit the GRO J 1655-40 spectrum. 
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Fig. 4. — IC342 Source 1 ASCA GIS spectrum fitted with the Schwarzschild disk model 
(i = 0°; top) and with the slim disk model (bottom). Residual is more conspicuous in the 
slim disk model below ~ 1.5 keV. 
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IC342 Source 1 Best— fit models 



Fig. 5. — IC342 Source 1 best-fit standard accretion disk model (face-on Schwarzschild disk; 
black) and slim disk model (red). Models removed of interstellar absorption are shown in 
dotted lines. 
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Fig. 6. — Spectral variation of IC342 Source 1 observe in September 1993. The same datasets 
as in Mizuno, Kubota and Makishima (2001) are used. The upper-panel indicates the varia- 
tion of mass and mass accretion rates obtained by applying the standard Schwarzschild disk 
model (GRAD model). Contours with dotted lines are from fitting in 0.5-10 keV, and those 
with solid lines are from 2-10 keV. The two contour levels indicate 68 % (1 a) and 90 % 
error regions for two parameters. The lower panel is obtained through application of the slim 
disk model in 2 - 10 keV. The region marked with gray indicates that the observed spectral 
variation is achieved with constant M and variable mass accretion rates. 





dM/dt (10 18 g s _1 ) 


Fig. 7. — Energy conversion efficiency of the slim disk used in the present paper and in 
Watarai, Mizuno and Mineshige (2001). The vertical broken lines indicate the critical mass 
accretion rates Me which gives the Eddington luminosity in the standard disk model (see 
appendix A). 
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